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We derive accurate formulas for thermal fluctuations in multilayer interferometric coating taking 
into account light propagation inside the coating. In particular, we calculate the reflected wave 
phase as a function of small displacements of the boundaries between the layers using transmission 
line model for interferometric coating and derive formula for spectral density of reflected phase in 
accordance with Fluctuation-Dissipation Theorem. We apply the developed approach for calculation 
of the spectral density of coating Brownian noise. 
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I. INTRODUCTION 

Thermal fluctuations in the mirror arc becoming signif- 
icant noise sources in the second generation gravitational- 
wave antennae (Advanced LIGO, Advanced VIRGO, HF 
GEO, TAMA) [1]. The pioneering articles on this issue 
dealt with Brownian fluctuations in the body of mirror 
[2-4]. Later the importance of thermoelastic noise [5] 
was realized especially for mirrors manufactured from 
sapphire. The physical reasons behind thermoelastic 
noise are provided by fundamental thermodynamic fluc- 
tuations of temperature manifesting through thermal ex- 
pansion. The same reason produces the thermorefractive 
noise [6] through the mechanism of relation between the 
refractive index and temperature. These results were ob- 
tained for model of infinite test mass (the mirror was 
considered occupying semi- infinite elastic space). Later 
these results were generalized for the finite-size mirror 
[4,7]. 

Soon the importance of thermal noise in the mirror 
coating was realized as the parameters of coating may dif- 
fer considerably from the mirror parameters bulk. Ther- 
moelastic noise in interferometric coating was calculated 
in [8, 9]. Later, following direction of H.J. Kimble [10] the 
potential of partial compensation of thermoelastic and 
thermorefractive noise in coatings was explored [11, 12]. 
However, today the Brownian noise in interferometric 
coating produces the main contribution into the noise 
spectrum [13-15] of gravitational wave antennae because 
the loss angles of substances used in coating (as usual, 
pair TaaOs, Si O2) are much greater than the loss angles 
in the same bulk materials. 

As a rule the calculation of Brownian and thermoelas- 
tic noise in coating is reduced to calculation of fluctu- 
ations of the total coating thickness. Strictly speaking, 
the light partially travels inside the coating, however, its 
power exponentially decreases with the depth when the 
light travels inside the coating. It allows considering the 
light to be perfectly reflected from the outer surface of 
the coating. Usually, the calculation of thermal noise 
in coating is provided under assumption that light is re- 
flected from the front surface of the coating without de- 
tailed analysis of the light propagating inside the coating 
[8, 9, 11-15]. 

The purpose of this article is to make accurate calcu- 
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FIG. 1: Top: Conventional mirror with interferometric coat- 
ing deposited on the face side. The coating consists of al- 
ternating quarter wavelength layers plus the outer half wave- 
length cap layer. Bottom: transmission line model with alter- 
nating quarter wavelength parts, pi = — P2 = — , Po = 1. 



lations of reflected wave phase taking into account prop- 
agation of light inside the coating. In particular, we cal- 
culate the reflected wave phase as a function of small dis- 
placements of the boundaries between the layers. In our 
analysis we use transmission line model for interferomet- 
ric coating [16]. In Sec. II we calculate the dependence 
of the phase of the wave reflected from the mirror as 
function of small displacements of boundaries between 
interferometric layers of coating. In Sec. HI we apply 
the formula obtained to calculate the spectral density of 
Brownian noise inside the coating to obtain numerical 
estimates. In Sec. IV we discuss obtained results. 



II. CALCULATION OF REFLECTED WAVE 
PHASE 

It is known that the interferometric coating consists 
of alternating quarter wave length layers with different 
refraction indices ni , n2 (ni > n2) as shown on Fig. 1. 
Speciflcally, below we assume that coating consists of 
TaaOs (ni) and fused silica (ria) layers, the outer layer 



(with positions zo, z^) is made from TQ2O5 and it is 
covered by additional half wave length layer (cap) which 
outer position is Zc- Such a coating is similar to trans- 
mission line consisting of alternating quarter wave length 
pieces with characteristic (wave) impedances pi = 1/n.i 
and P2 = l/fia [16]. 

We assume that the test mass (the body of the mir- 
ror) is manufactured from fused silica, hence, the wave 
impedance between zn and Zn+i is equal to pz — l/u-a- 
The mean position Zi of the surface of each layer (see 
Fig. 1) corresponds to the quarter wave length thick- 
ness of each layer, however, small position fluctuations 
Ci (caused by thermal fluctuations) produce fluctuations 
in reflected wave phase. In this section we calculate the 
reflected wave phase as a function of small displacements 
Ci using the successive approximation technique. 



A. Main formulas and the zeroth approximation 

In our consideration we ignore the optical losses in ma- 
terials of coating. We start with considering the last layer 
assuming that all fluctuation displacements Ci are equal 
to zero. We assume that the wave propagating in the 
positive direction of axis z is described as ~ e^^'^'^^^^'^'. 
The equivalent impedance at position z = Zn+i is equal 
to Zn+1 = Po = 1 (vacuum). It is convenient to define 
the amplitude reflectivy R as 



U,efl = RUi, 
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Obviously, the amplitude reflectivity Rn + i at z = Zn+i 
is equal to: 
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Now we can calculate the effective impedance Zn at posi- 
tion Zn — it allows considering the piece of transmission 
line between zn and Zn+i as a single impedance Zn 
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Here 9n + i is the exponent describing phase advance of 
wave traveling between positions zn and Zn+i , ^^,2 are 
the wave vectors, A is the optical wavelength in vacuum. 
Here and below superscript '°' refers to the particu- 
lar case of the zeroth approximation when fluctuations 
Cn, Cn + 1 are absent. We consider the test mass as one 
fused silica layer ("slightly" thicker). For convenience 
we assume its thickness being a fold to the quarter wave 
length (to meet the condition of anti resonance) — it 
means 0n + i = — 1- We make this assumption to make 
final formula more compact, without it the final result 



does not change practically but formulas looks more ak- 
ward. 

In the similar manner we successively calculate reflec- 
tivity Rn, impedance Zn-i and so on: 
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The total number N of layers is assumed to be an odd 
number. Hence, we can calculate the impedance at point 
z = Zo in the zero approximation: 
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Accounting for the half wave length cap layer gives that 

the effective impedance Zc at position Zc is equal to 

7(0) _ 7(0) 

Reflectivity and transparency. In case of small 
impedance Zo "C po one may easily calculate the effective 



amplitude reflectivity Ro and transparency To = W 1 — Rq 
Ro = ^^--1+^, To^^ (2.7) 

Po + Zo Po po 



In the zeroth approximation we have: 
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B. Fluctuations of layer positions in coating 

Now we take into account displacements Cj in position 
of each layer. We are interested in only linear terms of 



expansion Zo over Cj • Hence, we can calculate each term 
separately, i.e. we can calculate a term, for example, 
proportional to a certain fluctuation Cm putting all other 
positions to be equal to zero: Ci^m = 0- 




FIG. 2: For calculation of the dependence on displacement 

Co. 

First we assume that the cap layer is absent (later we 
will include it into considerations). We start by taking 
into account a small displacement Co assuming that all 
other displacements are equal to zero: Cj^^o = 0. Hence, 
we can account for the unperturbed piece of transmission 



r(0] 



see 



line between zi . . . Zn + i as a single impedance Z 
Fig. 2. 

The displacement Co influences the phase of reflected 
wave in two ways: a) through variation of effective 
impedance which, in turn, changes the reflectivity and 
b) through direct displacement of front surface which re- 
flects the wave (see formula (2.15) below). 

a). We calculate the perturbed impedance Zq: 
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Now we can estimate the perturbed reflectivity Rq: 
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b). Let the amplitude of the incident wave to be A 
then the complex amplitude B of the reflected wave is 
provided as 

B = RoA e^^'^^" (2.15) 



Now expanding exponent in (2.15) in series over Co and 
substituting (2.13) we obtain 
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We see that contribution of displacement Co into the 
phase of reflected wave is depressed by a small factor 

-(Pl/P2)^. 

In the same manner we can calculate contribution of 
fluctuation displacement in each layer Ci into the phase 
of reflected wave. In addition we can take into account 
fluctuations of position Cc of additional half wavelength 
layer (cap) made from fused silica. Then we can calculate 
the phase of the reflected wave expressing it as a function 
of small displacements Ci with account of the cap — see 
detailed calculations in Appendix A 
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FIG. 3: The coefficients Cj from formula (2.17) for a coat- 
ing with a cap. Red circles correspond to odd subscript and 
blue diamonds - to even ones, ni = 2.035 {Ta205),n2 = 
1.45 (SiOi). The number of TaaOs layers is assumed to be 
Ni — 20, and the number of Si02 layers — N2 = 19 plus cap. 
Subscript _i on the plot corresponds to ec. 

The formula for phase (^ of reflected wave is convenient 
to write through refractive indices (pi ~ 1/tl1) p2 = 
Vn2, n-i >rL2): 
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We see that fluctuation displacements Ci , Cs)--- (odd 
numbers) provide the main contribution into sum (2.17) 
whereas the input of displacements Co Co> C2) C4 • • • 
(even numbers) is negligible. The plot of coefficients ei 
is presented in Fig. 3. 

Note that formulas for coefficients e,, obtained in the 
cases with and without half wavelength cap arc practi- 
cally the same — compare formulas (A12) and plots in 
Fig. 8 in Appendix A with formulas (2.17) and plot in 
Fig. 3 correspondingly. 



III. CALCULATION OF COATING BROWNIAN 
NOISE 

We can apply the exact formulas (2.17) to calculate 
Brownian noise to coating. In accordance with Fluctu- 
ation Dissipation Theorem (FDT) [3, 17, 18], in order 
to calculate the spectral density Sx(cu) of fluctuations of 
variable X = ^ eiCi at the frequency cu, we have to ap- 
ply force Y.X €.1^0^^^^ acting at frequency cu so that the 
force eiFo is to be applied to position zi. Then we have 
to calculate the total dissipated power W and to find the 
spectral density using formula 



Sx(a)) = 



SIcbTW 



(3.1) 



In so doing we just calculate the spectral density of fluc- 
tuating variable X as there is no need to calculate addi- 
tional correlations between thickness fluctuations in dif- 
ferent layers. 

In contrast to previous approximate approaches [8, 9, 
11-14] the formulas (2.17) provide the option of direct 
explicit calculation of the thermal noise in coating. Recall 
that in above articles the following approximate formula 
was used for reflected wave phase 

•Papprox ^ •^^'iCc (.'J-2J 

In this case we have to apply force Foe^^^ to position z^ 
in order to calculate spectral density of variable Cc • 



Obviously, the use of explicit formula (2.17) will give 
smaller value of dissipated power relative to result of ap- 
proximate approach (3.2). Hence, the spectral density 
should be also smaller. 



A. Mirror as an infinite half space 

We compare the explicit and approximate calculations 
of Brownian noise of coating using the model of the mir- 
ror as a semi-inflnitc half space using approach developed 
for structural losses in [13, 14]. In accordance with FDT 
[3, 17, 18] we have to apply pressure eipo to the i-th 
boundary between layers: 
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and so on. We are interested in elastic energy Ui stored 
in each layer. 

First we write down formulas for strains u^j and stress 
Czz in outer layer of the mirror body (substrate) [4, 13, 
19] 
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Here A, \i are the Lame coefficients of substrate, which 
arc known to be expressed through the Young modulus 
Y and the Poisson ratio v as following 
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We assume that the number i = N corresponds to the 
surface of substrate (as shown in Fig. 1).. 

Now we calculate strains and stresses in coating layers. 
As usual we assume that the tangent strains in the i- 
th layer of coating (the positions of its boundaries are 
Zt_i , Zi) are equal to corresponding strains in substrate: 
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The last components u^^ of the normal strain in the i-th 
layer can be found from the known formula binding stress 



and strain tensors [19]: 



(3.8) 



The normal stress cr^^ • can be easily calculated as fol- 
lowing 
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FIG. 4: The coefficient Vi proportional to energy stored in 
the i-th layer from formula (2.17) for coating with a cap. Red 
circles correspond to odd subscript (riodd = 2.035 (Ta205)) 
and blue diamonds (ricvcn = 1-45 (Si02)) - to even ones. 
Number of TQ2O5 layers is assumed to be Ni = 20, and the 
number of Si02 layers — N2 = 19. Coefficent Vo on the plot 
corresponds to the cap layer. 
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Here (|)i is the loss angle of structural losses in the i- 
th layer, A denotes the optical wavelength in vacuum, 
thickness di of i-th layer is expressed through its refrac- 
tive index tli as di — A/4ni for each layer except cap, the 
multiplier [1 + 6oi] is introduced to account for the fact 
that the cap thickness is two times larger (as its width 
is equal to halj wavelength) . The plot of dimensionless 
coefficients Vi is presented on Fig. 4, the used numerical 
parameters arc presented in Table I. 



TABLE I: Parameters used for numerical calculations. 



Parameter 


substrate 


Ta205 layer 


Si02 layer 


T, K 




290 




A, m 




1.064 X 10-^ 




NcM 


- 


20 


1 9 + cap 


n 


1.45 


2.035 


1.45 


Y, Pa 


72 X 10' 


140 X 10' 


72 X 10' 


V 


0.17 


0.23 


0.17 


4) 


4x 10"^° 


2 X lO-'* 


4 X 10"^ 



Now we may compare the estimate of accurate formula 
for spectral density Sx of Brownian noise with approxi- 
mation S^^ used previously [13, 14]. Numerical calcula- 
tions for parameters presented in Table I gives: 



Manipulating with (3.8, 3.9) we write down a useful 
formula: 

. 2{A + M.)Ii-AiI 

u^.AiU.., A,.A_^L^_^ (3.10) 

and to calculate the elastic energy Ui in the i-th layer 
with thickness di through its volume density Wi using 
the known formulas: 
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Now we calculate the spectral density of variable X 
(3.1) using the fact that for the structural losses the 
power dissipated in the i-th layer can be expressed as 
Wi = Ui4)icu: 
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0.0539 



(3.14) 



We see that accurate calculation gives slightly less value 
of spectral density. However, the difference is quite mod- 
est — about 5 percents only. We may qualitatively ex- 
plain it by the fact that the tagent strains Urr, 11.4,4, are 
the same for all layers and they make considerable con- 
tribution into elastic energy. The account of different 
coefficients ei change only strain u^^, which is smaller 
than tagent strains for the majority layers. 



B. Mirror as a finite cylinder 

We repeated our calculations of coating Brownian noise 
for mirror as a finite cylinder using results of [4, 7, 15]. 
For estimates we used the parameters listed in Table I. 
For fused silica cylinder planned in third generation 
laser gravitational antenna (Einstein Telescope, radius 
R — 0.30 m, height H = 0.31 m, w = 0.12 m) we have 
got 
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For cylindric test mass planned in Advanced LIGO (ra- 
dius R = 0.17 m, height H = 0.20 m, w = 0.06 m) we 
have got 



Scc-Sx 



Sx 



0.0569 



(3.16) 
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We see that the difference is practically the same as for 
the model of mirror as an infinite half space. 



effective tansmittanccs of corresponding coatings), to cal- 
culate the power dissipated in coating and, finally, to cal- 
culate spectral density Sapp(uj) using formula (3.1). We 
have also calculated spectral density Saccl^^) using our 
approach through coefficients ei. We have found that 
Sacc(tu) is smaller than Sapp(uj) by about 17%: 



0.172 
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C. Double mirror (Khalili etalon) 

The number of layers used in conventional coating is 
large enough (about 40) — it is the reason why explicit 
and approximate formulas give close numerical results as 
estimates (3.14, 3.15, 3.16) show. However, explicit for- 
mulas may give considerably different results in double 
mirror [20]. Recall, in conventional mirror the fiuctua- 
tions of thickness of each layer in the coating are trans- 
formed into phase fluctuations of rcficcted wave in two 
ways. First, each layer makes contribution into varia- 
tion of front position of coating (position of first layer) 
and each layer makes approximately equal contribution. 
(It is worth underlying that we are interested in posi- 
tion of front surface of mirror relative to its center of the 
mass.) Second, the fluctuations of the layer thickness 
vary the pathlength of light traveling inside it and only 
several first layers make the main contribution whereas 
inner layers make exponentially small contribution. The 
first effect is much bigger than the second and usually 
only it is taken into account. 

The idea of double mirror, put forward by F. Khalili 
[20] (now the term "Khalili etalon" is frequently used), 
is to displace the part of layers from the front surface to 
the rear surface of the mirror. In this case the thickness 
fluctuations of layers on the rear surface do not make 
contribution into fluctuations of front surface relative to 
mirror's center of the mass. And number of layers on the 
front surface may be smaller than total number of layers. 
Obviously, the explicit calculation of Brownian noise of 
the front layers have to give considerably less value of the 
spectral density as compared to approximate one due to 
relatively small number of layers. Obviously, there are 
precisely layers on front surface which make the main 
contribution into the coating Brownian noise of double 
mirror. 

As an example, we have applied approach presented 
in this paper for Khalili etalon in order to estimate dif- 
ference between two approaches. We assume that the 
coating on front surface consists of 3 layers of TaiOs and 
2 layers of SiOa plus cap (also manufactured from fused 
silica) and coating on rear surface — 17 layers of TQ2O5 
and 17 layers of SiOi. The other parameters were taken 
the same. One may calculate coating Brownian noise 
in Khalili etalon by two ways. In traditional (approx- 
imate) approach one has to apply corresponding forces 
£f Fo and £rFo forces to front and rear surfaces of cylin- 
dric mirror (the coefficients £f and £r are calculated from 



Here we used the parameters from Table I and Advanced 
LIGO (radius R = 0.17 m, height H = 0.20 m, w = 
0.06 m). 



IV. CONCLUSION 

The formulas (2.17) may be applied to the explicit 
calculation of Brownian, thcrmoclastic, thermorefractive 
noise in the coating. These formulas are especially im- 
portant for the thermal noise compensation first proposed 
by J. Kimble [10] and then dcmonstarted for thcrmoclas- 
tic and thermorefractive noises [11, 12]. Recall that in 
[11, 12] consideration was based on approximation (3.2) 
but we hope that formulas (2.17) will allow to formulate 
explicit recommendation for thickness of each layer. 

In order to do it we have to rewrite formula for re- 
flected wave phase 4^ in form useful for the calculation 
of thermorefractive noise. For this case we have to as- 
sume that layers positions are not fluctuated Cj = but 
there are fluctuations of light path lengths due to varia- 
tion of refractive index in each layer — see corresponding 
formulas (A14) in Appendix A. This formula does not 
account influence of cap because in accordance to recom- 
mendation of [12] the thickness of additional cap layer is 
a subject for optimization. 

Recently non-quarter wavelength coating was sug- 
gested by I. Pinto [21]. As loss angle in TaaOs is much 
bigger than in SiOa one can use, for example, 1/8 wave- 
length layer of TaiOs and 3/8 wavelength layer of SiOi 
in order to reduce thermal coating noise. We plan to 
apply our approach for detail analysis of such coatings. 

We hope that proposed approach will be useful for de- 
tailed analysis of coating noise. 
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FIG. 5: For calculation of dependence on displacement Ci 



FIG. 6: For calculation of dependence on displacement ^2- 



Appendix A: Calculations of reflected wave phase 

Here we present detailed derivation of formulas (2.17) 
starting with taking into account of displacement Ci as- 
suming that all others are zero: Cj=,^i =0. Hence, we may 
consider the unperturbed piece of transmission line be- 
tween Z2 . . . Zn + ] as a single impedance Z2 — see Fig. 5. 
Again, the displacement Ci influences in two ways: a) 
variation of optical pathlength between zo and Zi and b) 
between zi and Z2 . We consider each of them separately. 

a). Note that influence of Ci on pathlength between zo 
and Zl we have almost calculated above. Indeed, calcu- 
lating influence of Co we have taken into account variation 
of exponent 9i and expanded 6^ ~ —(1 -|-2iki (Ci — Co)) 
keeping only term ~ Co- Now we may account displace- 
ment Ci just making substitution Co -^ ^ Ci in (2.12) 



Zo =Z[,°' -ikCiPo (1 -a^) 



(Al) 



h). We calculate variation of Zi by Ci and then recal- 
culate it into Zq: 



p(0) _ Z2_-P2 _ 1 - tXN-1 _ Po /Pi 

zi^p^l^^, e^^-(i-2ik2Ci), 
I — K2O2 

P2 [^ i-l<2Cl(1 -Kn-I 



N-1 



Fig. 6. Again we note that displacement C2 influences in 
two ways: variation of optical pathlength a) between zi 
and Z2 and b) between Z2 and Z3. 

a). We may account displacement C2 just making sub- 
stitution Ci -^ — C2 in (A2): 

„2 

Zo=z|,°'-ikpoC2{1-4_i)^. (A4) 

P2 

b). Wc calculate variation of Z2 by C2 and then recal- 
culate it into Zq: 



cxn-1 



aN-i 



Zo = |i=z(°'+ik2Ci(1-a^N_i)^ 

Zl p2 

= Z^°'+ikpoCi{1-a^_,)^. 

P2 



(A2) 



Collecting (Al) and (A2) we obtain the total contribution 
of Ci in impedance Zo 



7 7(°' ^r, ivr i P^~P2 
Zo = Zq + Po ikCi ( 2 — 

P2 



(A3) 



Now we take into account only displacement C2 as- 
suming that all others are zero: Cj7^2 — 0. Hence, we 
may consider the unperturbed piece of transmission line 



0) 



between Z3...ZN+1 as a single impedance Z3 — see 



2 2 / \ N-3 

P2 _ PW Pi 



,(0) 



Po VP2 



> Po, 



-,(0) 1 / \ N-2 

Z3 -pi 1 -aN-2 Po / Pi \ 

Z' ' +pi 1 +aN-2 P2 VP2/ 



Z2 = P 



Z2^P 



1+R39^ 

1-^363 

aN-2 + (1 -aN-2)il<iC2 



I, 0i^-(1-2ikiC3), 



1 +(1 "aN-2)ikiC2 
Pi oin-2 + Po ikC2(1 - aN-2), 



(0) 



P2 



Z2 = Z^"' + PoikC2(1-a^_2) 



Pi 



(A5) 
(A6) 



Collecting (A4) and (A6) we obtain the total contribution 
of C2 in impedance Zo 

Zo = Zo°' + Po ikC2 (o^N-i - ^^-2) \~j] = 

.4..,,„.,,(|)(£i)--k_^,., 

By the same consideration we may take into account 
displacements C3 and C4: 



zo.zr+P.iU3(ji)^(Hi^,+ 



(A8) 



I 2 2^ / \ 2 / \ 2N— 4 

and write down the final formula for effective impedance 
of transmission line as a function of small displacements 



Ci: 



Zo = Po 



N + l 



ikCoPo 1 



Po Pi 

P2 VP2 



(AlOa) 



+ p,ikCi (^^) + 



Po ikCs 



+ Po i.l<C4 



Pi VPi/ VP2 

\ 2 / 2 2 

Pi \ / Pi - P2 



P2 



9l 



(AlOc) 



r 2 1\ / \ 2 / \ 2N-4 

(Pi -P2) (?0\ /Pi 



Pi 



P2 y V p2 



--^-'iii)""'(^ 



Po ikCam 



(Pf -P2) /Po\^ / Pi 



Pi 



P2 y V p2 



(AlOd) 
- (AlOc) 

2N-2m 



62 



63 



64 



e2m-l 



e2T 



^N + l 



Tlt-m 



m 



nt -m 






ni-n^. 



2 2 



2 /'^li 
ni 

X 2 

Hi 

2 /t^2 



2N-2 



2N-4 



-) 

2 /^lli 

ni 



2TR-2 



n; 



2 /^ 
ni 



2N-2m 



N-1 



(A12d) 

(A12e) 

(Al2f) 



(A12g) 



(A12h) 



We see that coefficients eim witii even index are rela- 
tively small and their contribution may be omitted. The 
main contribution is made by coefficients eim— 1 with odd 
subscripts — see plot on Fig. 8. 
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FIG. 7: For calculation of equivalent impedance of the coating 
with a half wavelength cap. 

The formula for the phase cf) of reflected wave complex 
amplitude B (analog of (2.15)) is more informative one: 

B = ^°~P° Ae^'""^" c^ 4°'a(1 + 14)), (Alia) 
Zo + Po 



N + l 

ict) = 2ik Y ej Cj- 

j=0 



(AUb) 



One may rewrite these formulas through refractive in- 
dices (pi = 1/ni, p2 = l/ni, TLi > ni): 



N + l 

tcj) = 2ik ^ ej Cj , 

j=0 

/ \ 2N 
2 / Tl2 

eo = n2 — 

VTl-1 



ei 



n^-n2 



(Al2a) 

(A12b) 
(A12c) 



FIG. 8: The coefficients £j from formulas (A12). Red circles 
correspond to odd subscript and blue diamonds - to even 
ones. U] = 2.035 (Ta205),n2 = 1.45 (SiOa). Number of 
TQ2O5 layers is assumed to be Ni — 20, and number of SiOz 
layers — N2 = 19. 

We may rewrite formulas (A12) for the reflected wave 
phase 4) in the form useful for the calculation of ther- 
morefr active noise. For this case we have to assume that 
layers positions are not fluctuated Cj = but there are 
fluctuations of path lengths due to variation of refractive 
index in each layer: 

e? =exp[2ik(ni,2 + 6Tii,2)(zj+i -Zj)] - 



[1 + TtAj] , Aj 



(A13) 



Here Aj is the relative fluctuation of refraction index in 
j-th layer. Finally the phase of reflected wave has the fol- 
lowing dependence on fluctuations of refractive indexes: 



N + l 



14) = l7t Y_ £j ^j ) ^i = — - 
j=i 



_ Suj 
n 



(Al4a) 
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FIG. 9: The coefficients Cj from formuias (A14). ni = 
2.035 (TQ205),n2 = 1.45 (Si02). Number of TQ2O5 layers 
is assumed to be Ni = 20, and number of St02 layers — 
N2 = 19. 
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1 friz 



n^ \ni 



^-ni 



2 ("[^ 
ni 



2(N-1 



2(N-2 



1 / \ 3 / / \ 2(N-3) 



£n 



1 f riz 



TLi \n^ 



N-l 



1 -n^ 



2 1"^ 

n-i 





(A14g) 



We present plot £j in Fig. 9. 

Account of a cap. To account for the outer half wave- 
length layer (cap) we have to consider the transmission 
line model shown in Fig. 7. In the zeroth approximation 



we put Cc = Co = and obtain Z 



(0) 



In the first 



order approximation we use the expansion in series for 



00 to account for layer positions fluctuations: 

e^~i+2ik2(Co-Cc), 



P2 



2Zo + (Z|,"'-p2)2ilc2(Co-Cc 



2p2 



.(0) 



p2)2ik2{Co-Cc 



= Zo +ik2(£,o-£,c 




(0)i2 



P2 



P2 



(A15) 



(A16) 



Now using (AlO) and ki^2 — l^Po/Pi,2 we rewrite: 
+ PoikCc I 1 



(Al7a) 



Po / Pl_ 

P2 VP2 



N + 1 



-o«o^(f)(i;^~ 



Po 1.1<Cl 



Pi -P2 
Pi 



2 2 / 2 \ / \ 2N— 2 

p.ikC2^i^(%)(^ 
P2 VP2/ VP2 



^ 2 / 2 2 

•T / I Pi \ /^Pl ^P2 

+ Po ikCs — ' ' 

.P2 



Pi 



+ (A17b) 

f (A17c) 

(A17d) 
+ (Al7e) 
(A17f) 



f 1 2^/ \2/ x2N— 4 



Pi 



, 2rri— 2 / 7 2 



P2 



Pi 



Po ikC2m 



(Pf -P2) (Po\^ I Pi 



Pi 



P2 / V P2 



(A17g) 
- (A17h) 

2N-2m 



From these formulas one can obtain expressions (2.17) to 
expand into series the reflected wave B 



B = 



^c + Po 
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